Universality and the zero-bias conductance of the single-electron transistor 
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The thermal dependence of the electrical conductance of the single-electron transistor (SET) in the zero-bias 
Kondo regime is discussed. An exact mapping to the universal curve for the symmetric Anderson model is 
established. Linear, the mapping is parametrized by the Kondo temperature and the charge in the Kondo cloud. 
Illustrative numerical renormalization-group results, in excellent agreement with the mapping, are presented. 



Nearly five decades ago, Anderson conceived a Hamilto- 
nian to describe the interaction between a magnetic impu- 
rity and otherwise free conduction electrons.— Once a daunt- 
ing theoretical challenge, the Anderson Hamiltonian yielded 
to an essentially exact numerical diagonalizationf2 followed 
by an exact analytical diagonalization.^''* From these and 
alternative approaches, physical properties were extracted, 
which eased the interpretation of experimental dataj^ the- 
oretical results provided unifying views of apparently un- 
related phenomena;^ and quantitative comparisons brought 
forth novel perceptions.— 

The last ten years were especially fruitful. Parallel ad- 
vances in scanning tunneling spectroscopy and in the fab- 
rication of nanostructured semiconductor devices enhanced 
the interest in transport properties. In 
both areas, numerous experimental breakthroughs and the- 
oretical analyses were reported, and the Anderson Hamil- 
tonian proved spectacularly successful in more than one 
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occasionisli^ 

Notwithstanding the substantial volume of exact results, 
certain aspects of the model remain obscure. Consider uni- 
versality, a concept important in its own right and by virtue 
of its diverse applications. Universal relations serve as bench- 
marks checking the accuracy of numerical data; as resources 
promoting the convergence of theoretical findings; and as in- 
struments bridging the gap between the theorist's tablet and 
the laboratory logbook. The conditions under which the 
Anderson model exhibits universal thermodynamical prop- 
erties were identifiediSiM Although one expects all proper- 
ties of the model to be universal in the same domain, few 
firm results for the dynamical and transport properties can be 
found in librariesisa Costi's et al's early effort showed that the 
transport coefficients for the symmetric Anderson model are 
universalis^ For asymmetric models — even ones that display 
universal thermodynamical properties — , nonetheless, the uni- 
versal curves fail to fit the numerical data, the disagreement 
growing with the (particle-hole) asymmetry. 

Puzzled by such contrasts, we have conducted a systematic 
study of the transport properties for the Anderson Hamilto- 
nian. We combined analytical and numerical-renormalization 
group (NRG) tools and paid special attention to universal- 
ity. In a preliminary report^^ we have discussed an Anderson 
model for a quantum dot side-coupled to a quantum wire, a de- 



vice comprising two conduction paths whose transport prop- 
erties are marked by interference i2^i2^i22i22ia2 Notwithstanding 
the constructive or destructive effects, we have been able to 
identify universal behavior throughout the Kondo regime, the 
parametrical domain favoring the formation of a magnetic mo- 
ment at the quantum dot and its progressive screening by the 
conduction electrons as the temperature is lowered past the 
scale set by the Kondo temperature Tk- Specifically, we found 
the thermal dependence of the conductance to map linearly 
onto a universal function of the temperature T scaled by the 
Kondo temperature Tk- The mapping is itself universal, i. e., it 
depends on a single physical property, the ground-state phase 
shift 5, into which the contributions from all model parameters 
are lumped. 

This report examines the alternative experimental set-up in 
which a quantum dot or molecule, instead of side-coupled 
to, is embedded in the conduction ^a.i\\ Pr^-^^r'^r'~^^r''^r'^ We 
show that the thermal dependence of the conductance maps 
onto the same universal function. Although linear, the map- 
ping now depends explicitly on a model parameter — an exter- 
nal potential applied to the conduction electrons — and hence 
contrasts with the conclusion in our previous report. This 
dependence accounts for distinctions between the transport 
properties in the embedded and side-coupled arrangements. 
At high temperatures, for instance, potentials appropriately 
applied to the conduction electrons in the side-coupled geom- 
etry drive the conductance from low values up to the ballistic 
limit 02 - 2e~/h. If the quantum dot is embedded in the con- 
duction path, by contrast, the high-temperature conductance 
is pinned at low values and virtually insensitive to potentials 
applied to the conduction electrons. Our analysis shows that, 
in the embedded configuration, the screening charge in the 
Kondo cloud parametrizes the mapping to the universal con- 
ductance curve. Since that charge is always close to unity, the 
mapping is never far from the identity, with maximum relative 
deviations around 20%. 

Our presentation focuses the mapping between the SET and 
the universal conductances. As illustrations we will present 
the results of a few Numerical Renormalization Group (NRG) 
runs. A discussion of the numerics, a comprehensive survey of 
the Kondo regime, and the comparison with the side-coupled 
geometry will be deferred to another report. 

The text is divided in five Sections, more technical aspects 
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of the analysis having been confined to the three Appendices. 
Section H] defines the model. Section HI] derives an expres- 
sion relating the conductance to the spectral density of the 
quantum dot level. Section |III| is dedicated to universality, 
and Section IIVI to the fixed points of the model Hamilto- 
nian and to an extension of Langreth's exact expression for 
the ground-state spectral density. Section |V] then shows that, 
in the Kondo regime, the thermal dependence of the conduc- 
tion can be mapped onto the symmetric-SET universal con- 
ductance. Finally, Section|Vl]collects our conclusions. 



I. SINGLE-ELECTRON TRANSISTOR 

Figure[T]depicts a single-electron transistor (SET), the pro- 
totypical example of embedding. The subject of numerous ex- 
perimental studies, the SET comprises two independent con- 
duction bands coupled by a localized level. 




FIG. 1: Single-electron transistor. A quantum dot (circle) bridges 
two quantum wires (rectangles). A gate potential controls the dot 
energy, while the symmetric potentials V„, shift the the energy of the 
wire orbitals close to the dot. 

Qualitatively, the physics of Fig. [T| was understood long 
before the first device was developed. Figure |2] displays the 
spectrum of the SET Hamiltonian H for zero coupling. The 
dot levels being then decoupled from the conduction bands, 
the eigenstates and eigenvalues of H can be labeled by the dot 
quantum numbers. For simplicity, we will only refer to the dot 
occupation n^. For fixed n^, the product of the lowest dot state 
by the conduction-band ground state is shown as a bold dash. 
The gray levels above it represent the excited states consistent 
with the same rid label. 

A small transition amplitude V between the quantum dot 
and the wires is sufficient to modify this picture. The ampli- 
tude y couples strongly each gray level to the degenerate or 
nearly degenerate states in the neighboring columns. Excep- 
tions are the lowest levels in the column labeled «j - N m 
Fig. |2] which are energetically distant from their neighbors 
and thus remain unperturbed to first order in the coupling. 
At low temperatures, with small in comparison with the 
energy 6& separating the ground state from the closest level 
in the neighboring columns, the dot occupation is frozen at 
rid - N, a. constraint that raises the Coulomb blockade against 
conduction through the dot. 

Adjustment of the gate potential Vd in Fig. [T] brings down 
the blockade. The potential shifts the energies of the dot 
levels and can be tuned to the condition 6& ^ 0, to im- 
pose degeneracy between the bold dashes in the iid - N and 
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FIG. 2: SET energies in the weak-coupling limit. The dot-level oc- 
cupation nj labels the energies. For each nj, the bold dash represents 
the conduction-band ground state, while the thinner lines represent 
excitations. The coupling between the dot and the two quantum wires 
mixes each level to the neighboring columns. 

lid - N + I columns. An infinitesimal bias is then sufficient 
to induce electronic flow between the wires throught the dot. 
The conductance peak, we see, at gate potentials Vd such that 
the ground-state expectation value of iid is half-integer, e. g., 
m rid |Q> = + 1/2 as 5£ in Fig.|2] 

Each peak identifies a resonance at the Fermi level. As the 
gate voltage is swept past SS - 0, the ground-state occupation 
changes rapidly from iid - N io rid — N + \, and as required 
by the Friedel sum rule, so does the ground-state phase shift. 
At moderately low temperatures, for thermal energies smaller 
than the average spacing between the bold dashes in the figure, 
the plot of the conductance as a function of the gate voltage is 
a succession of peaks. Data collected in the laboratory at mod- 
erately low temperatures do display a sequence of resonances. 
At very low temperatures, however, the pattern changes to a 
sequence of intervals alternating between insulation and con- 
duction. 

The conducting plateaus are due to the Kondo effect. For 
gate voltages corresponding to odd ground-state dot occupa- 
tions, the magnetic moment of the resulting dot spin inter- 
acts antiferromagnetically with the conduction electrons. As 
the device is cooled past the Kondo temperature, the screen- 
ing of the moment creates the Kondo resonance, a spiked en- 
hancement of the density of states pinned at the Fermi level. 
Notwithstanding the Coulomb blockade, the pinned resonance 
allows conduction. 



II. ANDERSON MODEL 

A variant of the Anderson Hamiltonian encapsulates the 
physics of the device in Fig. [T] A spin degenerate level Cd 
represents the dot level, and two structureless half-filled con- 
duction bands, labeled L (left) and R (right), represent the 
two quantum wires. The L (R) wire comprises state 
(Qr) with energies defined by the linear dispersion relation 
eit - {k - kf)vf (0 < k < 2kf), so that the bandwidth is 
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2D — Ivfkp. The per-particle, per-spin density of conduction 
states is p = 1 1 2D, and we will let A = D/N denote the energy 
splittings in the conduction bands. The model Hamiltonian is 
then the sum of three terms, H = //„ + Hj + //,„/, where the 
first term describes the wires: 



ka kc/a 



(1) 



with an intra-wire scattering potential W, fixed by the poten- 
tial y„. in Fig.[Tl and a - L,R. The Hamiltonian Hd describes 
the dot: 



Hd = BdUd + UridfUdi, 



(2) 



with a dot energy e,/ controlled by the gate potential Vd in 
Fig- HI and the Hamiltonian H^d couples the wires to the dot: 



H. 



V 



y[2N 



^(cLq + H.c.). 

ka 



(3) 



where /(e) is the Fermi function; 



1 + ji^p^W^ 



(9) 



is the width T - npV^ of the Cd level, here renormalized by 
the scattering potential W; and 

Pd(e, T)^j--rYj ^l<«l 4 - e) (10) 

is the spectral density for the dot level. Here |m) and |«) are 
eigenstates of Ha with eigenvalues £■„, and £„, respectively, 
e,„n = E,„ - E„, and X is the partition function for the Hamil- 
tonian Ha ■ 

As one would expect, given that the odd Hamiltonian Hb 
commutes with Cd, only the eigenvalues and eigenvectors of 
Ha are needed to compute the right-hand side of Eqs. ([8]l and 
( fTOl ). The following discussion will hence focus the even 
Hamiltonian, Eq. (|5]), which is equivalent to the conventional 
spin-degenerate Anderson Hamiltonianji 



A. Parity 

To exploit the inversion symmetry of Fig. [l] we define the 
normalized even (aj.) and odd (b/i) operators 



ak = -^(ckL + Cku); 

bk - —=(CkL - CkR). 

V2 



(4a) 
(4b) 



The projection of the model Hamiltonian on the basis of the 
a^'s and bk^ sphts it in two decoupled pieces, H - Ha + Hb, 
where 

Ha^Yj ^'^-^l^'^- + ^^ofo + y(fo^d + H. c.) + Hd, (5) 

A- 

where we have introduced the traditional NRG shorthand 

fo = J]ak/^/N, (6) 



and 



HB^J]ekblbk+^Yj^lb. 



N 



(7) 



kq 



B. Conductance 

The odd Hamiltonian Hb is decoupled from the quantum 
dot. It is, moreover, quadratic, and hence easily diagonaliz- 
able. Appendix |C] determines its spectrum, analyzes the re- 
sponse of the conduction and dot electrons to the application 
of an infinitesimal bias and turns the result into the following 
Linear Response expression for the conductance: 



G(r) = ^2^rv 



f 

J-D 



Pd(e, T) 



dm 



de 



de. 



(8) 



C. Characteristic energies 

Four characteristic energies govern the physical properties 
of the Anderson Hamiltonian. Two of them are displayed in 
Fig. [3j the energy -e^ needed to remove an electron from 
the dot level; and the energy ed + U needed to add an elec- 
tron to the level. The particle-hole transformation q c|; 

Cd — > -c^ swaps the two energies, so that, the transformed dot 
Hamiltonian is given by the right-hand side of Eq. (|2]i with 
ed -(ed + U). 

If 2ed + U - 0, the dot Hamiltonian remains invariant un- 
der the particle-hole transformation. If, in addition, W - 0, 
Eq. ^ reduces to the symmetric Hamiltonian 

hI^Yj ^i<4''i< + ^(fo^d + H. c.) - y (n^t - ndif. (11) 
k ^ 

With y 0, two other energies arise: the level width Fw 
[Eq. ©I and the Kondo energy kBTfc, given by 



Tk ~ ^JpJexp(-l/pJ), 



(12) 



where / is the antiferromagnetic interaction between the con- 
duction electrons and the dot magnetic moment^^ 



p7 = 2 



U 



7T\Sd\ Sd + U 



(13) 



In the Kondo regime, thermal and excitation energies are 
much smaller than min(|ef/|, e,; + U). In Fig. [3] only the low- 
est levels in the central columns are energetically accessible. 
The energy Tw, associated with transitions from the central 
to the external columns in the figure (i. e., with and 

— > transitions) becomes inoperant. Instead, at very low 
excitation and thermal energies, smaller than the Kondo en- 
ergy kBTx, the dot spin binds antiferromagnetically to the con- 
duction spins. In Fig. [3] the lowest states in the left and right 
central columns hybridize to constitute a Kondo singlet. 



4 



■a 



1' 



w 









— -t- 1 


\ -4- 



ed + U 



LM 



FL 



FIG. 3: Spectrum of the spin-degenerate Anderson model, displayed 
as in Fig. [2] In the weak-coupling limit, the eigenstates are labeled 
by the occupation nj and spin component of the dot configuration 
displayed at the bottom. For V Q, each level in the left and 
right columns hybridizes with nearly degenerate levels in the cen- 
tral columns and acquires the width Yw in Eq. {9}. At low energies, 
the levels in the two central columns combine into a singlet and ac- 
quire a width Tf^ ~ /rsT^.The vertical arrows near the right border 
mark the domains of the LM and FL fixed points. 



III. UNIVERSALITY 
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FIG. 4: Thermal dependences of the conductance for two sets 
of model paramters, obtained from Eqs. l[8) and l IlQI l. The solid 
line depicts the universal conductance curve^^ for the symmetric 
Hamiltonian Jilt . Here, it was computed with F = 0. 1 Z) and 
U = 3 D. The temperatures were scaled by the Kondo tempera- 
ture Tk = 2.4 X 10"* D/kg, fixed by the requirement G(Tk) = O.Sffi- 
The dashed curve is the conductance for the Hamiltonian Q with 
F = OAD, U = 3D, £j = -0.3 D, and W = 0, which yielded 
Tk = 4.0 X 10"' D. To keep the data within the temperature range 
kgT < OAD, the dashed plot stops at T = 25 Tx- The horizontal 
arrows pointing to the vertical axes indicate the corresponding fixed- 
point conductances, given by Eqs. l |22at and l |22b| l. 



The concepts recapitulated in Section III CI emerged over 
three decades ago, with the first accurate computation of the 
magnetic susceptibility of the Anderson model,- long be- 
fore the first essentially exact computation of the conduc- 
tance. A particularly important result in Costi's, Hewson's, 
and Zlatic's survey of transport propertiesS^ is the the thermal 
dependence of the conductance for the symmetric Hamilto- 
nian H^, the universal curve G^{T/Tk), depicted by the solid 
line in Fig. |4l For ksT «c D and any pair (F, U) satisfying 
F <K t/ in Eq. ( fTTT ). proper adjustment of the Kondo temper- 
ature Tk gives a conductance curve G{T/Tk) that reproduces 

In Fig. m for instance, the solid line was computed from 
the eigenvalues and eigenvectors of with F = OA D and 
U ^ 3D. The definition G(Tk) = 0.5^2 yielded the Kondo 
temperature Tk - 2. Ax 10"^ D. When the calculation was 
repeated for U - 0.6 D and the same F, the Kondo temper- 
ature grew four orders of magnitude, to Tk - 2.2 x 10"- D. 
Still, for ksT < 0.1 D, the plot of 0(7 /Tk) resulted indistin- 
guishable from the solid curve. While Tk is model-parameter 
dependent, G{T/Tk) is not. 

Particle-hole asymmetry drives G away from G^ . For 
U + 2€d it or W ?i 0, the universal curve G^{T/Tk) no 
longer matches G(T/Tk). An example is the dashed curve in 
Fig.H calculated with F = 0.1 D, f/ = 3 D, e,/ = -0.3 D, 
and W = 0. The definition G(Tk) = 0.5^2, which in this 
case yields T^: = 4 x 10"^ D, forces the solid and the dashed 
lines to agree at T - Tk', the conductance for the asymmet- 
ric model nonetheless undershoots (overshoots) the universal 
curve for T < Tk {T > Tk). To reconcile this discrepancy 



with the concept of universality, the following sections rely 
on renormalization-group concepts. 

IV. FIXED POINTS 

Renormalization-goup theory probes the spectrum of 
Hamiltonians in search of characteristic energies and scaling 
invariances. The wire Hamiltonian ([1]), for instance, exhibits 
a single, trivial characteristic energy: the conduction band- 
width 2D. For energies e <sc D, therefore, its spectrum is 
invariant under the scaling transformation — > AH„, for ar- 
bitrary scaling parameter A > 1 . Accordingly, for e <sc D, the 
wire Hamiltonian is a stable fixed point of the renormalization- 
group transformation in Ref . 

Latent in the Anderson Hamiltonian (|5]), by contrast, are 
the four nontrivial characteristic energies discussed in Sec- 
tion lll CI Part of the spectrum of Ha lies close to fixed points; 
the remainder is in transition ranges. In the vicinity of a fixed 
point, the spectrum remains approximately invariant under 
scaling; in the transition intervals, the eigenvalues are compa- 
rable to one or more characteristic energies and hence change 
rapidly under scale transformations. In particular, the portion 
of the spectrum pertinent to the Kondo regime comprises two 
lines of fixed points and a crossover region. 

For given thermal or excitation energy £, the inequality 
max(fi, Fiy) <«; min(|e^|, + U, D) defines the Kondo regime. 
As Fig. [3] shows, the dot occupation is then nearly unitary. In 
the energy range ksTK «: £ «; minde^l, + U,D), which 
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is removed from characteristic energies, the Hamihonian Ha 
is near the Local Moment fixed point (LM). At very low en- 
ergies, & <s: ksTK, i. e., below the energy scale defined by 
the narrow set of levels at the center of Fig. [3] the spectrum 
becomes asymptotically invariant under scaling as the Hamil- 
tonian approaches the Frozen Level fixed-point (FL). In the 
intermediate region S ~ kgTK, the Hamiltonian crosses over 
from the LM to the FL. 



A. Fixed-point Hamiltonians 

As the two central columns in Fig. [3] indicate, the LM is 
an unstable fixed-point consistent of a conduction band and a 
free spin-1/2 variable. In the FL, a singlet replaces the spin, 
and the Hamiltonian is equivalent to a conduction band — a 
stable fixed point. In their most general form, the fixed-point 
conduction bands mimic the wire Hamiltonian, i. e.. 



The FL eigenvalues are likewise uniformly phase-shifted, 



and 



(14) 



(15) 



with scattering potentials Wfi and Wlm dependent on V, W, 
U and e^. Equations ( fT4b and ( fTsT i identify two lines of fixed 
points, parametrized by Wlm and Wfl, respectively. 

The Schrieffer- Wolff transformation offers an approxima- 
tion for the LM potential: 



pWlm = pW + 2 



Tw 2Ed + U 
n\Ed\ Sd + U ' 



(16) 



For most applications, this expression is insufficiently accu- 
rate, and an NRG computation is necessary to determine Wlm 
and Wfl- The exception is the Hamiltonian (fTTT i. for which 
Wlm - 0, as required by particle-hole symmetry. 



f^kg,,gk. 



(20) 



where et — et- {6/7t)A. From the Friedel sum rule, it follows 
thalS 



S - 5lm - -■ 



For Ha - H^a' ™ particular, 5 - njl. 



(21) 



C. Conductance at the fixed points 

The LM is the fixed point to which the Anderson Hamilto- 
nian would come if F - 0. ForO < r min(|erf|,e^/ + [/, D), al- 
though the renormalization-group flow never reaches the LM, 
it brings Ha close to the fixed point. The substantial portion of 
the spectrum of Ha marked by the thin double-headed arrow 
in Fig.[3]is approximately described by the many -body eigen- 
values of H^^, and in the pertinent energy range, the physical 
properties of Ha and H*^^^ are approximately the same. Like- 
wise, at low temperatures, the properties of Ha approach those 

of//;^. 

The renormalization-group evolution of the Hamiltonian 
can be traced in the termal dependence of the conductance. 
As the temperature is reduced from T » Tk to T <s: Tk, each 
curve in Fig. |4] crosses over from a lower plateau to a higher 
one. The extension of Langreth's expression^! derived in Ap- 
pendix |B] determines the plateau conductances: 

Glm = Si sin^(5LM - 5w) = Gi cos-((J - 6w); (22a) 
Gfl = 02^m^{d-6w), (22b) 

where 6w is the ground-state phase shift for V = 0. According 
to the analysis in AppendixlAl 



B. Fixed-point phase shifts 

Appendix |A] diagonalizes the quadratic Hamiltoni- 
ans (O and ( fTSl ). For the LM, the diagonal form reads 



HIm = Yj ^f^slst^ 



with phase-shifted energies 



5lm . 

se = €{ A. 

7T 



(17) 



(18) 



At the LM, all conduction states are uniformly phase-shifted, 
with 



tan^LM = -ttpWlm- 



(19) 



For Ha - H^, in particular, 6lm - 0, and the low-energy 
eigenvalues e^. coincide with the e^. 



tan^iv — -npW. 



(23) 



The solid curve in Fig |4] was computed for Ha - H^a^ 
that 6w - 0, while the ground-state (i. e., FL) phase shift is 
5 - njl. According to Eqs. ( l22l i. Glm - and Gfl - Qi. in 
agreement with the plot. The ground-state phase shift for the 
dashed curve, extracted from the low-energy eigenvalues in 
the NRG run that generated it, is somewhat lower: 6 - 0.437r. 
Again b\N - 0, and the two horizontal arrows pointing to 
the vertical axes in Fig.|4]indicate the conductances predicted 
by Eqs. ( |22] |. Given the relatively high Kondo temperature 
{kBTK = 4 X 10^3 D) in this run, the condition ksT <s; £> re- 
stricts the curve to the range T < 25 Tk, so that, even at the 
highest temperature shown. Ha is relatively distant from the 
LM, and Eq. ( I22ab cannot be accurately checked. At low tem- 
peratures, however, the renormalization-group flow bringing 
Ha asymptotically close to H^^, the agreement with Eq. ( |22bl ) 
is excellent. 
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V. CROSSOVER 



and sum the result over q, to find that 



In the Kondo regime, the SchriefFer- Wolff transformation^^ 
brings the Anderson Hamiltonian Ha to the Kondo form 

k liv 

with J defined in Eq. ( fTsT l. 

To eliminate the scattering potential on the right-hand side, 
it is convenient to project Hj upon the basis of the eigenoper- 
ators gk of the LM, which yields^ 

= ^ sig\gc + 4>l^o-^y(pQy ■ S, (25) 



where Jw - J cos^ 6lm, and 



(26) 



In the symmetric case 6lm vanishes, and the operator 0o re- 
duces to /o. 

The second term on the right-hand side of Eq. ( |25T l drives 
the Hamiltonian from the LM to the PL. Along the result- 
ing trajectory, the eigenvalues of Hj scale with 
Let Tk and Tk < Tk be the Kondo temperatures corre- 
spondig to two sets of model parameters in the Kondo regime: 
M = {F, W, U, Bd) and M = {f , W, U, Sd], to which corre- 
spond the antiferromagnetic couplings J and /, respectively. 
If |m) is an eigenvector of Hj with eigenvalue £„,, then a cor- 
responding eigenvector |m) of Hj, the scaling image of \m), 
can always be found, with the same quantum numbers and 
eigenvalue E,„ such that E„,/Tk = E,„/Tk- 

The matrix elements of any linear combination of the oper- 
ators gk are moreover universal. Given two eigenstates \m} and 
|n) of Hj and their scaling images \m} and \h}, then the matrix 
elements of <po, for example, are equal: {m\ (po \n} - {m\ 0o !«)■ 
Likewise, the matrix elements of the operator 

are universal: {m\ (pi \n) - {m\ (pi |«). 



A. Thermal dependence of the conductance 

By contrast, the matrix elements {m\ Cd \n) on the right-hand 
side of Eq. ( fTOb are non-universal. Even at the lowest energies, 
as Eq. ( IB 12b shows, they depends explicitly on the model pa- 
rameters. To discuss universal properties, therefore, we must 
relate them to universal matrix elements, such as {m\ (po \n), 
{m\ (pi \n}, or {m\g( |«). As a first step towards that goal, we 
evaluate the commutator 



[HA,al] 



V 



W 

'n 



(28) 



lHA,4]-^f!+Vcl + w4. 



Here we have defined another shorthand 



-5 V 9 



(29) 



(30) 



Equation ( |29l ) relates the matrix elements of between two 
(low energy) eigenstates \m} and |«) of Ha to those of the op- 
erators /() and fi : 

V{m\ cl \n) = (£,„ - E„ - W){m\ ^ |«) - V3D<m| // |«). (31) 

In the Kondo regime, with max(£'„,, £„) <s: D, the first two 
terms within the parentheses on the right-hand side can be 
dropped. 

In the symmetric case, since /o (fi) coincides with 0o {(pi), 
Eq. (|3TI ) shows that the product V{m\ Cd |n) is universal, in line 
with the firmly established notion that YpdielksTK, T/T^), 
and G^{T/Tk) are universal functions To discuss asym- 
metric Hamiltonians, we have to relate the operators /o and 
fi to (p{) and (pi. This is done in Appendix I A 21 which shows 
that, in the Kondo regime, a linear transformation with model- 
parameter dependent coefficients relates the matrix elements 
of both /o and fi to those of (po and (pi. When Eq. ( IA18I ) is 
substituted for fo and fi on the right-side of Eq. OTT i. it results 
that 



^jnp^w{m\ |«) = Q'o<'«l <Po !«> + ai{m\(p\\n}. (32) 

Here, the constants ao and ai are combinations of the (un- 
known) linear coefficients on the right-hand side of Eq. ( IA18I ), 
the parameter W on the right-hand side of Eq. OTI ). and the ra- 
tio -sJnprw/V, by which we multiplied Eq. (ISTT i to shorten the 
following algebra. 

Substitution in Eq. ( fTOl i yields an expression relating the 
spectral density pd to universal functions: 

npTwPdie, T) = Q'oPo(e, T) + a\pi(e, T) + aoaiP(oi)(e, T), 

(33) 

where 



Pj 



-\{n\(Pj\m)\' 



zm' 

x6(E,„-E„-e) (; = 0,1), (34) 



and 



P(01) 



(e,r) =2 



(<m| 4>l \n){n\ (pi \m) + c. c.) 



zm 

x6{E„,-E„-e). 



(35) 



Next, we substitute Eq. ( 1331 ) on the right-hand side of 
Eq. to split the conduction into three pieces: 

G{T) = alGoiT) + ajdiT) + (^oaiG(oi)(r), (36) 
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where 



P J-D 



dfie) 



de 



de 0- = 0,1), (37) 



and 



&2 r" 

G(or){T) = — P(Oi)(e, T) 

P J-D 



dm 



de 



de. (38) 



B. Universal contributions to the conductance 

Given the universaUty of the energies E,„ and of the ma- 
trix elements {m\ (pj |«) {j - 1 , 2) on the right-hand sides of 
Eqs. (|34] | and (1351 ), we see that the spectral densities p/(e, T) 
(j - 0,1), and P(oi)(e, r) are universal. Inspection of the 
right-hand sides of Eqs. ( |37||38t shows that the functions Gj 
(j = 0, 1) and G(oi) are likewise universal. To compute them, 
we are free to consider any convenient Kondo-regime Hamil- 
tonian. 

Particle-hole symmetry makes especially convenient. 
To show that the cross terms make no contribution to the con- 
ductance, i. e., that G(oi)(r) = 0, we only have to notice that, 
while leaving unchanged, the particle-hole transformation 



8k sl ( 



1. e., 



a,) reverses the sign of the 



product of matrix elements {m\ \n){n\ (pj \m) + c. c. on the 
right-hand side of Eq. ( [35] l. We see that p^^^^(e, T) is an odd 
function of e, so that the integral on the right-hand side of 
Eqs ( |38] | vanishes. 

To evaluate Go and Gi, we start out from the closed form 
resulting from the diagrammatic expansion (in the coupling 
V) of the conduction-electron retarded Green's function for 
the symmetric Hamiltonian: 



tm - 



f\e)6,, + ^Gf (e)G^(e)G[°\ 



(39) 



where is the retarded dot-level Green's function for the 



symmetric Hamiltonian, and 



Gf\e) 



1 



e-ek + ir] 



(40) 



is the free conduction-electron retarded Green's function. 

From G^j,, it is a simple matter to obtain the spectral densi- 
ties on the right-hand side of Eq. ( |33] ): 



Po(e, T) 



(41) 



and 



Pi(e, T) 



ttND^ 



3 J]e,ei-Gf,,(e). (42) 



kk' 



To compute the conductances at temperatures T satisfying 
ksT «: D, we only need the spectral densities for e «; D. 



It is appropriate, therefore, to expand the right-hand side of 
Eq. ( |40b to linear order in e/D: 



l'f\e)^^-i7T6(e-et) 



(e « D). 



(43) 



The sums over momenta on the right-hand side of 
Eqs. (|4TI ) and (l42l i are then easily computed. Among the 
resulting terms, only the even powers of e contribute to the 
integral on the right-hand side of Eq. (l37b . To compute the 
conductance to Olik^T I D)~] we hence neglect the terms of 
0{e/D). Equation ( HTT i then gives 



po(e, T)^p-Kprp',(e, T), 



(44) 



equivalent to an expression obtained in Ref.— . 

Substitution of this result for po on the right-hand side of 
Eq. ( |37] i establishes a simple relation between the universal 
function Go and the universal conduction for the symmetric 
Hamiltonian: 



GQ(T)^g2'G'(T). 



(45) 



Equation ( |45] l becomes exact, asymptotically, at low tem- 
peratures. The deviations, of 0[{kBT / D)^], are insifignicant. 
As an illustration, the open circles in Fig.|5]show NRG data for 
the conductance Go{T), Eq. dJTl i. in excellent agreement with 
the solid line representing the right-hand side of Eq. 
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FIG. 5: NRG results for the thermal dependence of the auxiliary con- 
ductance Go(T), associated with the spectral density for the operator 
00. The open circles show Eq. ( I37t for j = 0, computed for the sym- 
metric Hamiltonian with the displayed model parameters. The solid 
line is the right-hand side of Eq. ( 145b . i. e., the universal curve in 
Fig.|4]subtracted from the quantum conductance ^2- 

To the same accuracy, we can neglect the 0{e/D) terms re- 
sulting from the summation on the right-hand side of Eq. (l42l) . 
which yields 



pi(e,r)= ^p^(e,r). 
n 



(46) 



Equation dJTl i then shows that Gi is also related to the con- 
ductance for the symmetric Hamiltonian: 



Gi(T) 



(47) 



C. Mapping to the universal conductance 

The combination of Eqs. (l45T l and ( l47b with the result 
G(oi) = reduces Eq. ( |36] l to the equality 



G(T)^al{g2-G'{T)) 



y\-pG'(T) 
n 



(48) 



To determine the coefficients ao and ai, we need only com- 
pare the right-hand side with the fixed-point expressions for 
the conductance. At the LM, - 0, and Eq. (I22al i shows 
that al = cos^((5 - 6w)- At the FL, = §2, and Eq. ( l22b] i 
shows that (6/7r^)Q'j = sin~(5 - Sw)- These two results turn 
Eq. (|48] ) into the mapping 



(49) 



D. Illustrative numerical results 

Equation ( fT2b off'ers an approximation for Tj^, and 
Eqs. ( fTSI l. (fT9] l and (l2Tl i provide an approximation for the 
ground-state phase shift 6. These estimates are far from the 
accuracy needed to fit numerical or experimental data. In the 
laboratory, Tk and 6 - 6w are adjustable parameters; the for- 
mer, in particular, is determined by the condition G{T£) - 
Q^I2^^2iliiL In the computer office, the two unknown param- 
eters on the right-hand side of Eq. ( |49l ) can can be extracted 
from the conductance itself, or from other properties of the 
model Hamiltonian. The phase shift 6 is most easily obtained 
from the ground-state eigenvalues of Hf,. To determine the 
Kondo temperature Tk, it has been traditional to fit the ther- 
mal dependence of the magnetic susceptibility x(T^) with the 
universal curve for kgTxiT ITk)^ Here, however, we pre- 
fer the laboratory definition, which insures that both sides of 
Eq. (Ell) vanish at T = Tk- 

Figure|6]displays the results of two NRG runs for the same 
assymetric parameters U - 3 D, e - -0.3 D, F = 0.1 D, with 
two scattering potentials W - Q and W - -0.6 D. The open 
circles reproduce the dashed curve in Fig.lD The particle-hole 
asymmetry, combined with the relatively high ratio between 
the dot width and excitation energy, F/|e| = 1/3, place the 
model Hamiltonian close to the border of the Kondo domain. 
Although W = 0, the ground-state phase shift deviates sig- 
nificantly from n/2: from the FL eigenvalues generated by 
the NRG diagonalization of the model Hamiltonian, we find 
6 = 0.437r. The solid curve through the center of the circles is 
a plot of Eq. ^ with 6-6w = OASn and keTK^^x 10'^ D. 

The scattering potential W - -0.6 D reduces the Kondo 
temperature and raises the FL conductance G{T = 0). The 
former shrinks to keTK = 2.2 x 10"** D, while the latter rises 
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FIG. 6: Numerical data for the temperature dependence of the con- 
ductance, compared to Eq. l |49t . The open circles and triangles show 
the NRG computed conductances for the indicated model parame- 
ters. The solid lines represent the mapping , with ground-phase shifts 
calculated from the FL eigenvalues of the model Hamiltonian and Tk 
determined by the condition G(Tk) = ^2/2. The small disagreement 
between the triangles and the solid line above kgT = 10"" D is due to 
the relatively large irrelevant operators introduced by the scattering 
potential W, whose contribution to G decays in proportion to kgT ID. 



to nearly §2- Both changes are due to the reduced dot width 
Yw = cos 5wT. Here Eq. © yields pTw = 0.74pr = 0.074, 
and the resulting smaller antiferromagnetic coupling ( fT3b 
brings the Kondo temperature (fT2t down five orders of mag- 
nitude. 

The diminished Kondo temperature indicates that the scat- 
tering potential has pushed the model Hamiltonian deeper in- 
side the Kondo regime. Other indications are the minute high- 
temperature conductance; the nearly ballistic low-temperature 
conductance; and the overall similarity between G{T/Tk) and 
the solid line in Fig.|4] 



£. Discussion 



In the Kondo regime, Eqs. ( I22al ) and i22H fix the high- 
and the low-temperature conductances, respectively. Equa- 
tion ( |49l ) shows that the universal function G^{T/Tk) con- 
trols the monotonic transition between the two limits. For 
W = 0, in particular, the fixed-point values depend only on 
the ground-state phase shift 6 and are symmetric with respect 
to §2/2: Glm - ^2 cos^ 5 and Gfl = ^2 sin" 5. Thus, de- 
pending on 6, the transition from Gim to G fi can be steeper 
or flatter. Since 6 can never depart much from :7r/2 in the 
Kondo regime, the argument of the trigonometric function 
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on the right-hand side of Eq. ( |49] l can never depart substan- 
tially from n, and as indicated by the two curves in Fig. IH 
G(T/Tk) ~ G^{T/Tk) + 20%. By contrast with this crude 
estimate, the mapping ( |49l ) gives excellent agreement with the 
circles in Fig.|6] 

The wire potential W narrows the dot level and displaces 
the ground-state phase shift. Depending on the sign and mag- 
nitude of W, the phase shift can take any value in its domain 
of definition -n/2 < 6 < n/l. In the Kondo regime, the 
Friedel sum rule nonetheless prevents the difference 6 - 6w 
from straying away from n/l. All effects considered, the scat- 
tering potential W displaces the conductance curve towards 
the symmetric limit G(T/Tk) = G^T/Tk). 

These findings are in line with the experimentally estab- 
lished notion that, in the Kondo regime, SET conductances 
always decay with temperature , ^■^■^O-^^ A brief comparison be- 
tween this behavior and that of the side-coupled devic o^^'^^ 
seems appropriate. As demonstrated in Ref. a linear map- 
ping analogous to Eq. ( |49] l can be established between the 
side-coupled conductance and G^{T ITk)', in that case, how- 
ever, the coefficient relating the two functions is indepen- 
dent of 6w and hence free from the constraint imposed by the 
Friedel sum rule. Under a sufficiently strong wire potential, its 
sign can be reversed. Thus, the thermal dependence of Gsc is 
tunable:^'' a wire potential can turn a monotonically increasing 
function into a monotonically decreasing one. The embedded 
geometry of Fig.[T]is much less sensitive to W. 

The parameter 6w is (tt times) the charge induced under the 
wire electrodes by the potential W. According to the Friedel 
sum rule,^^ the difference 5 - i5iy is the charge of the Kondo 
cloud, the additional charge that piles up at the wire tips sur- 
rounding the dot as the temperature is lowered past Tk- Neu- 
trality makes the charge of the Kondo cloud equal to the dot 
occupancy. Since the symmetric condition «</ = 1 maximizes 
the the low-temperature conductance, one expects G{T - 0) 
to be ballistic for 2(<5 - 6w) = tt, a conclusion in agreement 
with Eq. (l49b . Since the screening charge is always nearly uni- 
tary, one expects the low-temperature conductance to be close 
to the conductance quantum, in agreement with the plots in 
Fig. in 



ence 6 - 5^/. Both are quantities of physical significance. Ac- 
cording to the Friedel sum rule, the phase shift difference is 
(7r/2 times) the screening charge surrounding the dot at low 
temperatures. 

In summary, we have derived an exact expression relating 
the SET conductance in the Kondo regime to the universal 
conductance function for the symmetric Anderson Hamilto- 
nian. A subsequent report will exploit this result in an attempt 
to offer a unified view of an NRG survey of conductance in 
the Kondo regime. 
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APPENDIX A: PROPERTIES OF THE FIXED-POINT 
HAMILTONIANS 

1. Diagonalization 

The LM and FL are described by conduction-band Hamil- 
tonians of the form 



k 

We want to bring H* to the diagonal form 



where 



(Al) 



(A2) 



(A3) 



To this end, we compare the expressions for the commutator 
[§f , H*\ obtained from Eqs. ( lAll) and (IA2I) . from which it fol- 
lows that 



VI. CONCLUSIONS 



acq 



1 W 

Be - eq N 



z 



atk- 



(A4) 



Our central result, Eq. ( |49] | maps the conductance in the 
embedded geometry onto the universal conductance for the 
symmetric Anderson model (fTTI ). Different from the univer- 
sal result for the side-coupled geometry, the mapping depends 
explicitly on the potential W applied to the wire. Section IVl 
showed that, in the Kondo regime, the Friedel sum rule an- 
chors the the argument of the cosine on the right-hand side of 
Eq. ( |49] ) to the vicinity of tt; it results that G{T /Tk) reproduces 
semiquantitative^ the universal function G^{T/Tx)- 

At the quantitative level, Eq. ( |49] l affords comparison with 
experimental data collected anywhere in the Kondo regime. 
For that purpose, its linearity is particularly convenient. Once 
fitted to a set of experimental points, the mapping determines 
the Kondo temperature Tk, as well as the phase shift differ- 



Summation of both sides over q then leads to the eigenvalue 
condition: 



W* 1 



N 



N ^E[- eg 



(A5) 



Inspection of this equality shows that, with exception of a 
split-off energy, which makes 0{l/N) contributions to the 
low-energy properties, the sc are shifted by less than A from 
the Q. We therefore refer to the closest conduction energy ee 
to label each eigenvalue and define its phase shift 6( with the 
expression 



£{ = e{ 6{. 

n 



(A6) 
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This definition substituted in Eq. ( IA5I ), a Sommerfeld- 
Watson transformation^ determines the sum on the right- 
hand side: 



-7 — 



C 1 

-npcotdf+p f 



■ de. (A7) 



Substitution in Eq. ( IA5I ) results in an expression for the phase 
shifts: 



TrpW* 7T J_o EC - 



■ de. 



(A8) 



At low energies, the contribution of the last term on the right- 
hand side, of 0{e/D), can be neglected, which shows that the 
phase shift becomes uniform: 



tan 5 = -npW* . 



(A9) 



To determine the coefficients acq, we square both sides of 
Eq. (IA4b and sum the result over q: 



(AlO) 



To evaluate the sum on the left-hand side, we differentiate 
Eq. iAli with respect to £[, which yields, with relative error 
0(1 /N), 



A?2 Zj _ e^^)2 [sin 6c I 



(All) 



The sum on the left-hand side of Eq. ( lAlOl ) being unitary, 
Eq. dAllI ) shows that 



W* a[k sin 6t 



(A12) 



the negative phase insuring that ffu — » 1 for W* 0. Equa- 
tion iA4i then gives 

A sindr 

ae, = [St « D). (A13) 



2. Energy moments of the matrix elements of the 
eigenoperators ge 

The Hamiltonian dAll ) diagonalized, we turn our attention 
to the following energy moments 

Mi^* = -^^(^y\mUd«> 0^ = 0,1,...). (A14) 

Chiefly important are M,5,^* s {m\ <po \n} and M„',|' = 
{m\ 4>\ \n)l VS; the other moments, as shown below, are pro- 
portional to M„'J,' = {m\ (pi |«). Since the M//„, are universal, to 
evaluate them it is sufficient to consider the symmetric Hamil- 
tonian ( fTTT i. for which the phase shift 6im - 0, so that gi,, e/^, 
(f)o, and (pi coincide with at, et, /o, and /i, respectively. 



From Eq. ( fTTI ). we then have that 



V 



[ge,H^] = S(gc + —=Cd. 

yN 



(A15) 



With the shortand G,„„ = (E,„ - E„)/D, the multiplication of 
both sides by (sf/D)''^^ followed by summation over £ leads 
to the coupled recursive relations 



M, 



M, 



ip) 

mn 

(p) 



SmnM, 



(p-1) 



V , 



-{m\ Cd |n> 



(P = 1,3,...); 
(P = 2,4,...). 



Reduced to a matrix equation, this system is easily solved. 
The result is 



Mln^&L - -{m\ Cd \n) 1 1 + 2^ rij^ 



p 
y_. 
p 



P-1 



M,J£',„„ <m| Cd |«> > 

n i—i r 



(p = odd) 
(p - even) 



where the primed sum is restricted to odd r's. 

The pertinent energies satisfy the condition kgT <s: D, 
which implies S,„„ 1. It is therefore safe to discard the 
terms proportional to &,„„ and its powers. Within this approx- 
imation, the only nonzero even moment is M,*",', and all odd 
moments are proportional to {m\ Cd |«). It follows that all the 
odd moments are proportional to M, 

{m\ 4>i |«) 



mn = {m\ (pi |«): 



M, 



(p) 







(p = odd) 
(P = 2,4,. 



(A16) 



.) 



An orthonormal basis describing the conduction band (IA2I) 
can be constructed from the definition 



^V^^ZW^^ (P-0,h...), 



(A17) 



where Pp{e) denotes a Legendre polynomial. 

According to Eq. ( IA16b . {m\ (pp |n) ~ {m\ (p\ \n) (p = 
3,5,...), while {m\ 4>p |«) = (p = 2, 4, . . .). This shows that 
the matrix element of any conduction operator is a linear com- 
bination of {m\ (po |«) and {m\ (pi \n). In particular 



(m| fi \n) = ^ aij{m\ cpj |«) (/ = 0, 1), 



(A18) 



where fii and fi are the operators defined by Eqs. ^ and dSOb , 
respectively, and the o-y (/, j - 0,1) are constants that depend 
on the model parameters. 



APPENDIX B: FIXED-POINT CONDUCTANCES 

This appendix derives an expression for the spectral density 
Pdie, T), defined by Eq. ( fTOl l, at the fixed points. The proce- 
dure is analogous to the one in AppendixlAl 
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From Eq. ( I28l l we obtain an expression for the matrix ele- 
ment of Uq between two low-energy eigenstates \m) and \n) of 
eigenstates of Ha '■ 



{m\ flj \n) 



1 



V 



■{m\c^\n) 



W 



1 



■{m\Y^al\n). (Bl) 



Summation of both sides over q leads to an expression for the 
matrix element in the last term on the right-hand side: 

{m\ 2 a], \n) (1 - WS,„„) = ^/NV{m\ 4 \n)S,„„, (B2) 



where 



which brings Eq. (IB II) to the form 

{m\ c[ \n) 



{m\ al \n) = 



V 



^/N(E„, - E„ - e,) 1 - WS„ 



(B3) 



(B4) 



Consider now this equality at one of the two fixed points, 
LM or FL. The fixed-point Hamiltonian has then the quadratic 
form ( IA2I) . which defines the complete basis of the operators 
g(. The matrix element (m|^j, \n) vanishes unless \m) = g\\n), 
which implies E„, - E„ + Sf. At a fixed point, therefore, 
the sum on the right-hand side of Eq. ( IB3I ) reduces to that in 
Eq. dATli, i. e.. 



<Smn - -Tip cot 6t, 



(B5) 



where we have disconsidered the last term on the right-hand 
side of Eq. ( IA7I ) because at a fixed point the ratio E(/D 
0. This result suggests that we introduce the phase shift 6w, 
defined by 



W = 



tan 5m 
Tip 



to simplify Eq. ( IB4l l: 



, , -I- , . yi>n\cl,\n} sin 5, cos 

{m\a \n) — . 

yfN(E„, - E„ - eq) sin(5. - 5w) 

In analogy with Eq. ( IA3b we can, moreover, write 
gc = ymcd + ^ ycqaq. 



with normalized coefficients: 



7^0 + 2 



1. 



(B6) 



(B7) 



(B8) 



(B9) 



At each fixed point, therefore, once squared, Eq. ( IB4l i reads 



N{se 



y^lsu^5^\ ^^^^^ 

- eq)- \sm((5, - 6w)l 



We divide both sides by A^, sum them over q, and substitute 
Eq. (lAl II) for the resulting sum on the right-hand side, to find 
that 

Substitution in the second term on the left-hand side of 
Eq. ( IB9I ) now shows that, with error 0(1 IN): 



\{m\c[\n)\^ 



1 sin {5t - 5w) 

NV^ TT^p^ COS^ Sw ' 



(B12) 



A the fixed points, the matrix elements are constants, de- 
pendent only on the phase shift and scattering potential. The 
spectral density p^, as one would expect, becomes indepen- 
dent of the temperature: 



equivalent to 



TT^p- COS'^ 



sin (5» - 6w) 



jrFcos- 5w 

W - ^ recovers the celebrated expression^ 



Pd{e) = 



sin^ 5* 
ttF 



(B14) 



(B15) 



More generally, however, to obtain the fixed-point spectral 
densities, we set 6* = (5 at the FL, and 6* - 6 - njl at the 
LM, from which it results that 



LM ^ C0S^(5 -6w) _ 

jtT cos^ 6w 
PL ^ sin^(5 - 6w) 
ttFcos^^w 



(B16a) 
(B16b) 



Substitution of Eqs. ( IB16al ) and ( IB16bb for pj on the right- 
hand side of Eq. dHJl leads to Eqs. (I22ab and ( I22bb . respectively. 



APPENDIX C: ZERO-BIAS CONDUCTANCE 

By contrast with the coupling to the impurity, which is in- 
dependent of the odd operators b); defined by Eq. ( l4b] i. the 
Hamiltonian describing a bias voltage couples to the bt's. Pre- 
liminary to the discussion of the conductance, it is therefore 
convenient to derive results for the b^^'s analogous to those in 
Appendix [A] Specifically, given the formal equivalence be- 
tween Eqs. ^ and ( lAll ). we can follow the steps in that Ap- 
pendix to write Hb in the diagonal form 



with 



k 



(CI) 



(C2) 
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and to derive a result analogous to Eq ( IA12l i. At low ener- 
gies, in particular, i. e., for \ei\ D, the eigenvalues ec are 
uniformly phase shifted by 6w [see Eq. ( |B6I )1. and 



= cos (5,1 



(C3) 



Multiplication of both sides by g[ and summation over E then 
shows that 

bk \h) = cos 5„ g( \h), (C4) 



for any pair |m), \h) of low-energy eigenstates of Hb- 

It is likewise convenient to compute the following commu- 
tator: 

[Ha, alh] = -^c>, + ^ Yj^4bk - albg), (C5) 

from which we see that, given two eigenstates |*P,„) and |«) of 
Ha with eigenvalues £„, and E„, respectively. 



W ^{'i'm\a\bk-alb,\^„) 



1. Current 



Em En 



(C6) 



To calculate the conductance, we can, for instance, examine 
the current flowing into the R wire; 



dqR 
dt 



(C7) 



1. e.. 



/ = -^[Ha, J] {aUk + blbk - (albk + H. c.))], (C8) 



which is equivalent to 

because summed over k, the last term on the right-hand side 
of Eq. iC5i vanishes. 



2. Conductance 

To induce a current, we add to the model Hamiltonian an in- 
finitesimal, slowly growing perturbation that lowers the chem- 
ical potential of the R wire relative to that of the L wire: 

//„ = Ati h^{t) = -e^ J] - 4c,^) e^'l\ (CIO) 



with an infinitesimal shift A/y. 

Projected on the basis of the a^s and bk^, reads 



and Eq. ( 



(Cll) 



shows that 



k>0 



{%„\clbk - b'lcM 



(C12) 

Standard Linear Response Theory relates /i^ to the conduc- 
tance: 



GiT) 



e-''^'"{%„\[I,h^it)]\%„)dt, (C13) 



where Z is the partition function at the temperature T. 

Comparison with Eq. dCllb shows that the operators at- 
within the parentheses on the right-hand side of that equality 
make no contribution to the conductance. We therefore define 



and rewrite Eq. (IC13 



GiT) 



e-''^'"{%M,hM)W,n}dt. (C15) 



Following the insertion of a completeness sum 2„ |«)(«| on 
the right-hand side of Eq. (IC15I ), straightforward manipula- 
tions lead to the familiar expression 



<^„,I4I%)<%|/V(0)I^,„) 
E,„ - E„ + ii] 



On the right-hand side, we now substitute Eq. ( IC12I ) 
[Eq. ( |CT4] |1 for li^ {h). This yields 

e2 y2 ^ (Qi',„\blcd\^„}{'i'n\clbk\'¥,„} 



m,n,k,q 



E„, - E„ + iri 
{'i'r,^clbk\'i'„}Qi'„\blca\^',„}] 



E-m Eifi + iff 
< . 

Em ~ En 



(C16) 



Aided by Eq. ( IC4I ). we can now trade the sum over the con- 
duction operators bk for a sum over the eigenoperators ge: 



G(T) = -i 



■ e Ely 
IhNpZ 



2 

n,n,e,e' 



{'VJglcd\%,){'i',,\cyge\^'m} 
E,„ - E„ + irj 

{'Vm\clgt\^„){'Vn\~g\,Cd\'V,n) 



E,„ - E„ + irj 



(CI 7a) 
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Since the g[ diagonalize Hb, only the terms with { - €' con- 
tribute to the sum on the right-hand side. We interchange the 
indices m and n in the second term within the parentheses on 
the right-hand side to show that 

G{T) = Yj e-^''"M'^m\clgt\'i',)\^6{E^ - E„). 

(C18) 

Since \^„,) - \m)\m), where \m) is an eigenstate of Ha 
(of the quadratic Hamiltonian Hb), the right-hand side splits 
into two coupled sums: 



G(T) 



~h NpZ 



2 e-P^-\{m\ Vc\ \n)\ME„, - E„ - q) 
:Y^e-l'^-{m\g(\nmg\\m). (C19) 



The second sum is equal to Zb[l - f(ep)], where /(e) is the 
Fermi function, and Zb, the partition function for the Hamil- 



tonian Hi,. The identity 

_ 1 ^/ 
"7(0 57 

then turns Eq. (IC19I ) into 
e~ ttFw 



(C20) 



G{T) 



hZa pN 



\{m\cl,\nt 



x6(E„,-E„-ee), (C21) 



where Xa is the partition function for the Hamiltonian H^- 

The definition ( fTOb of the spectral density pd{e, T) allows 
us to rewrite Eq. ( IC21l i as 

from which Eq. ^ follows. 
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